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COUP D'CEIL SUR L'ANALYSE DE FOURIER 



Embrasser I'analyse de Fourier d'un coup d'oeil est liors de ma portee. Je 
promenerai le regard sur I'liistoire et sur certains des termes en usage. L'his- 
toire est ancienne et elle est instructive. Je commencerai par Platon. 

Chez Platon, les mathematiques comprennent cinq parties : les nombres, 
les figures planes, les solides, la musique et I'astronomie. La musique et I'as- 
tronomie sont comme deux soeurs, liees par la musique des spheres. Et le 
programme astronomique de Platon est de rendre compte du mouvement des 
astres errants, les planetes, pour le rendre conforme a I'harmonie du monde et 
a la dignite des dieux. 

Le systeme de Ptolemee a realise ce programme, par une superposition de 
mouvements circulaires uniformes, le grand cycle et les epicycles. Et quoique 
le systeme de Ptolemee soit abandonne, la decomposition d'un mouvement en 
mouvements periodiques et I'utilisation des fonctions trigonometriques a ete 
une constante de I'astronomie jusqu'a nos jours. Aujourd'hui I'astrophysique 
double et depasse I'astronomie dans I'approche des phenomenes cosmiques. 
Mais la perle de I'astronomie qu'est la decouverte des exoplanetes repose sur 
I'analyse de Fourier des spectres emis par les etoiles. 



Texte des exposes faits aux Journees mathematiques X-UPS le 10 mai 2011. 
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Quant a la musique, son lien aux mathematiques est multiple. Le princi- 
pal est la decomposition d'un son en harmoniques et sa reconstitution comme 
somme d'harmoniques. C'est au sens propre I'analysc ct la synthcsc harmo- 
niques. A ccla s'attachcnt Ics gammcs, la distinction dcs timbres, Ics cordcs 
vibrantes et la grande controverse qui au 18° siecle a oppose Daniel Bernoulli 
a d'Alembert, Euler et Lagrange sur la possibilite de representer une fonction 
par une serie trigonomctrique. Aujourd'hui I'analyse par ondelettes complete 
ou supplante I'usage des fonctions trigonometriques pour tenir compte des 
echelles de temps aussi bien que des frequences, et le champ de I'analyse de 
Fourier s'etend par ses methodes comme par ses usages. 

Fourier, selon Riemann, est le premier a avoir compris completement la na- 
ture des series trigonometriques, en associant ce que j'appellerai I'analyse, les 
formules intcgrales donnant les coefRcients, et la synthese, la representation 
d'une fonction par une serie trigonometrique. C'est a cause de Riemann que 
nous parlous aujourd'hui de series de Fourier. En France, Fourier a ete long- 
temps meconnu. Arago, dans son eloge de Fourier, vante grandement le savant 
et le politique, explique I'importance de sa theorie analytique de la chaleur, 
mais ne dit pas un mot des scries trigonometriques, c'est-a-dire de I'outil que 
Fourier a forge pour calculer des solutions des equations integrales de la cha- 
leur. Fourier y accordait une grande importance et, contrairement a une idee 
regue, il a developpe cet outil en veritable theorie, exemples, applications et 
demonstrations, avec une grande rigueur. Mais il s'est heurte a I'incomprehen- 
sion persistante de Lagrange, et deux pages dans les manuscrits de Lagrange, 
que j'ai consultees et commentees, confirment que Fourier avait raison contre 
Lagrange. Mais Lagrange etait a I'epoque de Fourier le plus respecte des ma- 
thematiciens frangais, et son jugement negatif sur Fourier a traverse les siecles. 
Dans I'edition que je possede d'Encyclopaedia Universalis il n'y a pas d'article 
sur Fourier. 

Fourier attachait une portee universelle a ses formules. II precisait bien, et 
le premier, que la donnee d'une fonction etait celle de son domaine de defini- 
tion en meme temps que d'une loi ou une figure ; et, sans utiliser ces termes, il 
distinguait soigneusement les intervalles ouverts et les intervalles fernics. Mais, 
apres avoir multiplie les exemples et donne quelques preuves, il s'etait aventure 
a dire que toute fonction etait sujette a cette analyse et representable par une 
serie trigonometrique qui converge vers la fonction. Litteralement c'est faux. 
Pour appliquer les formules integrales, il faut que la fonction soit integrable ; et 
la convergence des series est un sujet difficile. Les premiers pas pour eclaircir 
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la question sont dus a Dirichlet, avec le premier theoreme general de conver- 
gence et le premier exemple de fonction non integrable dans la conception de 

I'cpoque. Tous les concepts d' integration, a commencer par I'intcgrale de Rie- 
mann, sont lies aux formules de Fourier et a leurs conditions de validitc. An 
dela de Riemann, on pense a Lebesgue, Denjoy, Laurent Schwartz. Quant a la 
convergence, elle est mise en question par les contre-exemples : une fonction 
continue dont la serie de Fourier diverge en un point (du Bois-Reymond) ou 
meme sur un ensemble donnc de mesure nulle (Kahane-Katznelson) , et une 
fonction integrable au sens de Lebesgue dont les serie de Fourier diverge par- 
tout (Kolmogorov) . Le theoreme de convergence de Carleson, inattendu a son 
epoque (1966), dit que la convergence a lieu presque partout quand la fonction 
est de carrc integrable ; on pent ameliorer cette condition, sans parvenir bien 
siir aux fonctions intcgrables. 

De meme que le concept d'integrale, c'est le concept de serie qui est en cause. 
Pourquoi s'attacher a la convergence, qui a d'ailleurs plusieurs significations 
des qu'on passe a plusieurs variables et a des series multiples, alors qu'il y a 
des procedes de sommation utilisables ? Le tournant est pris avec le theoreme 
de Fejer de 1900 : pour les fonctions continues, les moyennes arithmetiques 
des sommes partielles convergent. EUes convergent meme uniformement. La 
convergence dans les espaces fonctionnels apparait pen apres; les fonctions 
sont des points, les sommes partielles d'autres points, et I'espace des fonctions 
de carre integrable s'impose a I'attention avec la formule de Parseval, etablie 
dans ce cadre par Fatou, et surtout le theoreme de Riesz-Fischer, qui eta- 
blit I'isomorphisme isometrique de et de £^ par transformation de Fourier 
(1907). 

Fischer et Frederic Riesz, independamment, out pense a une geometrisation 
des espaces de fonctions, et le lemme fondamental pour la preuve de leur 
theoreme est le meme ; nous Tcxprimons aujourd'hui en disant que est 
complet. Mais cette formulation ne date que de Banach dans sa theorie des 
operations lineaires de 1930. II fallait jusque la une longue phrase pour le dire. 
C'est un exemple ou les definitions, tardives, viennent exprimer le sue d'une 
methode, avant de servir de base a de nouveaux developpements. 

Autre exemple, toujours tire de I'analyse de Fourier. Au lieu de et l^, 
Norbert Wiener s'est attache a et l^. Leurs transformees de Fourier sont 
un champ d 'etude toujours ouvert, avec des applications surprenantes en theo- 
rie du signal (Donoho, Candes etc.). Les problemes d'analyse et de synthese 
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y prennent un aspect different, plus algebrique : dans I'algebre des fonc- 
tions sommes de series trigonometriques absolument convergentes, les fonc- 

tions nuUes sur un ensemble donne forment un ideal ferme ; y en a t-il d'autres ? 
C'est bien la cas, coninie Malliavin I'a niontrc en 1959. Le point de depart est 
I'algebre de Wiener, qui se voit soit comme I'algebre multiplicative des fonc- 
tions sommes de series trigonometriques absolument convergentes, soit comme 
algebre de convolution La derniere expression est plus rapide, et la notion 
de convolution, si fondamentale en analyse, la rend tres parlante. Mais en 1930, 
et meme quand Laurent Schwartz a elabore sa theorie des distributions, on 
ne parlait pas encore de convolution ; c'etait Faltung en allemand, produit de 
composition en frangais. Dans le traite de Widder « Laplace transforms », qui 
date de 1941 , c'est comme « Stieltjes transforms » que sont introduites les 
convolutions de mesures, avec une note en bas de page qui indique le ternie de 
convolution comme une timide nouveautc parallclenient au terme de Faltung, 
aussi utilise en anglais. La convolution apparait partout, mais c'est Wiener qui 
a degage son caractere fondamental , et c'est pourquoi Faltung est devenu pour 
un temps I'expression de la notion. Aujourd'hui il est raisonnable de placer la 
convolution au depart d'un cours d'analyse de Fourier. 

S'agissant du vocabulaire, comment situer I'analyse harmonique? C'est un 
champ largement ouvert sur I'ensemble des mathematiques. Historiquement, 
on vient de voir ses relations avec les equations differentielles de la mecanique 
celeste, avec les equations aux derivees partielles des cordes vibrantes et de la 
chaleur, avec la theorie des fonctions d'une variable reelle, et on a evoque ses 
relations actuelles avec la statistique et le traitement des donnees. Le lien aux 
probabilites est ancien et profond. Le cadre de I'analyse harmonique commu- 
tative est celui des groupes abeliens localement compacts, et la dualite entre 
ces groupes est exprimee par la transformation de Fourier. La these de Tate 
a niontrc I'importance de cette approche en theorie des nonibres. L'analyse 
harmonique non commutative est liee aux groupes non commutatifs et a leurs 
representations. L'analyse harmonique abstraite part de la theorie de Gelfand 
des anneaux normes, autrement dits algcbrcs de Banach. L'analyse de Fou- 
rier classique traite de transformations intcgrales, transformation de Fourier 
d'abord, et aussi integrales singulieres. Ses usages dans toutes les sciences out 
ete multiphes par la transformation de Fourier rapide puis par les ondelettes 
et leurs variantes. II est clair que toute definition de l'analyse harmonique 
en serait une limitation injustifiee. Mais outre son etendue, on pent reperer 
des questions, des methodes, des theories qui, elles, peuvent etre identifiees et 
formalisees. 



ANALYSE ET SYNTHESE HARMONIQUES 



5 



Les formules de Fourier sont I'exemple de base. Leur enonce par Fourier 
exprime leur generalite, de fagon formellement incorrecte. EUes ne constituent 
pas un theoreme, mais elles ont engendre des theoremes et permis de defi- 
nir d'importantes notions. Mieux qu'un theoreme, elles ont constitue un pro- 
gramme, a savoir, donner des conditions de leur validite. Elles ont ensuite 
constitue un paradigme pour tous les developpements orthogonaux. Une rai- 
son de leur succes est sans doute qu 'elles etablissent un pont entre deux classes 
d'objets, en I'occurrence des fonctions et des suites. La dualite de Fourier est 
un modele de traduction d'un langage dans un autre, un trait commun a 
d'autres grands programmes. 

Le coup d'oeil pourrait se poursuivre, aussi bien sur I'histoire que sur I'actua- 
lite de I'analyse de Fourier. Ce sera pour une part I'objet des exposes suivants. 
En attendant, voici quelques references. 

Rene Spector (1970), dans I'article « Harmonique (analyse) » Encyclop<£dia 
Universalis, donne un bon apergu en quatre parties : series de Fourier, analyse 
et synthese harmonique, transformation de Fourier, groupes commutatifs 
localement compacts, avec un mot sur les groupes de Lie. 

Jean-Paul Pier (1990), dans L'analyse harmonique, son developpement his- 
torique (Masson) prend pour point de depart I'apparition de la theorie des 
groupes topologiques generaux, de la mesure de Haar et des algebres de Ba- 
nach, done il se limite pour I'essentiel au xx'' siecle, avec beaucoup de de- 
tails sur les questions qu'il traite, et une abondante bibliographie (40 pages) 
|Pie90] . 

Elias Stein et Rami Shakarchi, dans Fourier analysis, an introduction (2002) 
ont un point de vue tres different. Au depart, pas d'integrale de Lebesgue ni de 
groupe : integrale de Riemann, droite et cercle suffisent, puis M*^, en enfin Z^r, 
avec beaucoup d' applications |S-S03] . C'est le premier livre d'une serie qui 
se poursuit avec l'analyse complexe, la theorie de la mesure, I'integrale de 
Lebesgue et les espaces de Hilbert, puis les distributions de Schwartz et les 
probabilites, dans un foisonnement d 'applications. 

On voit la variete des approches. D'excellents livres ont paru dans les an- 
nees 1930. Je signale seulement le grand classique, le traite d'Antoni Zyg- 
mund Trigonometrical series qui, dans la seconde edition (1959) est devenu 
Trigonometric series ; la « troisieme » edition (2002) reproduit la seconde 
|Zyg35^ [ZSFp2t |Zyg59^ [ZS^ . 



Abstract harmonic analysis de L.Loomis (1953) a permis aux gens de ma 
generation de s'initier aux algebres de Banach (les anneaux normes d'Israel 
Gelfand). Le seul point commun aux livres de Zygmund et de Loomis est 
la theorie de Norbert Wiener des transformees de Fourier des fonctions inte- 
grables [Loo53| . 
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Sur ce sujet, qui a fait de grands progres au cours des annees 1950, Wal- 
ter Rudin fait le point avec Fourier transforms on groups (1962) jRud62| 
IRudQO] . Des progres ulterieurs sont exposes dans Essays in commutative har- 
monic analysis (1980) de Colin Graham et O.C.McGehee [ GM79] . 

Henry Helson (1983) litre simplement son livre Harmonic analysis et il se 
centre sur un petit nombre de questions, typiques et essentielles ; son but, 
dit-il, est la simplicite, et c'est un modele d'ecriture |Hel83|, IHellOj . 

Au regard de H. Helson, Tom Korner est foisonnant dans ses deux livres 
Fourier Analysis (1988) et Exercises (1993). Ces deux livres, savants et riches, 
pleins d'apergus et d'anecdotes, offrent une delectable promenade dans I'ana- 
lyse de Fourier |Kor88|, IKor93] . 

C'est aussi une promenade historique que Pierre-Gilles Lemarie-Rieusset et 
moi proposons dans Series de Fourier et ondelettes, en citant abondamment 
les grands ancetres, et en detaillant I'histoire recente des ondelettes [KLR98J. 

Je termine par le livre d'Yitzhak Katznelson, An introduction to harmonic 
analysis, parce que j'ai assiste a sa naissance dans le Languedoc et a son 
audience croissante aux Etats-Unis et dans le monde (1968-1976-2006). C'est 
beaucoup plus qu'une introduction, c'est I'expose clair et vivant de I'analyse 
harmonique par I'un de ses principaux acteurs |Kat68L IKat76L IKat04j . 



PARTIE I 
STATISTIQUE ET FOURIER, 
DANS L'HISTOIRE ET AUJOURD'HUI 



Ce sera une promenade historique, dans laquelle je pense m'attarder sur 
Fourier, mais qui aura comme point de depart la methode des moindres car- 
res et comme point d'arrivee I'echantillonnage parcimonieux, alias compressed 
sensing. 

Personnages 



Clairaut 1713-1765 
Lagrange 1736-1813 
Laplace 1749-1827 
Legendre 1736-1813 
Fourier 1768-1830 
Gauss 1777-1855 



Riemann 1826-1866 

E. Fischer 1875-1854 
Lebesgue 1876-1941 

F. Riesz 1880-1948 
N. Wiener 1894-1964 



Yves Meyer 
Stephane Mallat 
David Donoho 
Emmanuel Candes 
JPK 
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Le personnage principal sera Fourier, mais je parlerai aussi de ses aines, Le- 
gendre et Laplace, et de ses successeurs, parmi lesquels Riemann et Lebesgue. 
Dans la foison dcs theories et des rcsultats du 20'' siecle j'insisterai sur deux 
aspects : ce qui concerne i'^, I'espace des suites de carre somniable, qui a etc 
mis en evidence par Frederic Riesz et Ernst Fischer, et ce qui concerne 
I'espace des suites sommables, qui a ete decouvert par Norbert Wiener. La 
statistique sera I'occasion de voir i'^ et £^ s'introduire naturellement, avec les 
moindres carres de Legendre et compagnie, et le compressed sensing d'Emma- 
nuel Candes, David Donoho et compagnie. 

Pour ne pas risquer de manquer de temps en m'attardant en route, je com- 
mence par la fin du voyage : le theoreme de Candes, Romberg et Tao de 2006 
sur la rcconstitution d'un signal ou d'une image a partir de donnees appa- 
remment insuffisantes. En thcorie du signal, on enseigne que pour reconstituer 
une fonction dont la transformee de Fourier est portee par une bande de lar- 
geur L on a besoin d'un echantillonnage sur une progression arithmetique de 
raison a/L, a dependant de la normalisation. En limitant la fonction a un 
intervalle de longueur H, on a done besoin de HL/a points. Or void le theo- 
reme : 

On salt qu'un signal est constitue de T pics, sans connaitre leur position, 
et on a acces a sa transformee de Fourier. II s'agit de reconstituer le signal 
a partir de la transformee de Fourier de la maniere la plus economique. On 
identifie le signal a une fonction / portee par T points du groupe cyclique 
a A'^ points, N choisi assez grand. Ce groupe cyclique Z^v que je note G, est 
isomorphe a son dual G, le groupe des frequences, sur lequel est definie la 
transformee de Fourier /. On choisit au hasard, avec la probabilite naturelle, 
f2 frequences ui et on considere les f{oS). Supposons 

T<C^/\ogN avec C=1/23(M + 1). 

Alors on pent reconstituer / a partir des f{ijj) avec une probabilite 
ou c est une constante absolue. Et la reconstitution se fait par un procede 
d'analyse convexe, a savoir en minimisant la norme des fonctions definies 
sur G dont la transformee de Fourier coincide avec / sur les frequences 
choisies. On voit I'economie considerable : I'utilisation classique de la trans- 
formee de Fourier necessite A'^ points, et I'echantillonnage aleatoire le reduit a 
un multiple de TlogN. 

Rappelons que la norme £^ d'une fonction h definie sur G est la somme des 
valeurs absolues des h{g), g parcourant G. C'est bien une fonction convexe 
de h. 11 y a aujourd'hui des proccdcs varies pour traiter de tels problemes de 
minimisation, et je dirai un mot de leur histoire. 
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Remontons a 200 ans en arriere, plus precisement en 1805. Legendre a ecrit 
un memoire intitule « Nouvelles methodes pour la determination des orbites 

des cometes » ou il utilise un precede qu'il decrit soigncuscmcnt dans un ap- 
pcndicc dc quatrc pages intitule « Sur la mcthodc des moindrcs quarrcs » , ct il 
applique ensuite le meme precede a la mesure des degres du meridien terrestre, 
un sujet de grand interet a I'epoque. 

Get appcndicc est lumineux. En voici le debut : 

« Dans la plupart des questions ou il s'agit de tirer des mesures donnees par 
I'observateur les resultats les plus exacts qu'elles peuvent offrir, on est presque 
toujours conduit a un systeme d'equations de la forme 

E = a + bx + cy + fz + etc. 

dans lesquellcs a, b, c, f, etc. sont des coefficients connus, qui varicnt d'une 
equation a I'autre, et x,y,z, etc. sont des inconnucs qu'il faut determiner par 
la condition que la valeur de E se reduise, pour chaque equation, a une quantite 
ou nuUe ou tres petite. » 

Je resume la suite. En general, il y a plus d'equations que d'inconnues, 
et il n'est pas possible d'annulcr toutes les erreurs E,E',E" etc. Le precede 
consiste a rendre minimum la seninic des carres des erreurs. En derivant la 
somme des carres des erreurs par rapport a chacune des variables, on obtient 
des formes affines en nombre egal a celui des variables, et en les annulant on 
pent calculer ces variables, et aussi les erreurs E, E' , E" etc. Legendre precise 
qu'« on aura soin d'abreger tons les calculs, tant des multiplications que de 
la resolution, en n'admettant dans chaque operation que le nombre de chiffres 
entiers ou decimaux que pent exiger le degre d'approximation dont la question 
est susceptible. » 

Si parmi les erreurs ainsi calculees il en est d'aberrantes, on rejette les 
equations corrcspondantes et on reprend le calcul, dont la majeure partie (les 
multiplications) aura deja ete faite. 

Au prix de la multiplication des coefficients et de la resolution d'un systeme 
d'equations lineaires on a done un procede general et efficace pour realiser le 
programme indique par Legendre. Le procede prend une forme remarquable 
quand il n'y a qu'une variable x, qui represente un point sur la droite ou dans 
le plan ou I'espace, et des observations etc. Le minimum de la somme 

des carres distances au point x des points x',x",x"' etc. est atteint au contre 
de gravite de x',x",x"' etc. Comme le dit Legendre en conclusion, 

« on voit done que la methode des moindres carres fait connaitre, en quelque 
sortc, le centre autour duquel viennent se ranger tons les resultats fournis par 
I'experience, de maniere a s'en ecarter le moins qu'il est possible. » 
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II y eut une querelle de priorite entre Legendre et Gauss sur la methode 
des moindres carres. Gauss I'a publiee plus tard et sous certaines hypotheses 

sur la loi des erreurs (en gros, qu'cllcs soient gaussiennes) et il I'a utilisee sans 
la publicr avant 1805. Mais il scniblc bicn que Legendre, et d'autres, I'aient 
aussi utilisee sans la formaliser avant 1805. Stigler, dans son histoire de la 
statistique, consacre son premier chapitre a la methode des moindres carres 
exposee par Legendre (en donnant en illustration I'image d'un autre Legendre, 
cela est une tout autre histoire), avec un commentaire enthousiaste : 

« For stark clarity of exposition the presentation is unsurpassed ; it must be 
counted as one of the clearest and most elegant introductions to a new statis- 
tical method in the history of statistics. » 

Et surtout Stigler situe la methode parmi celles qui avaient cours prece- 

demment, a savoir les combinaisons ad hoc des equations pour se ramener a 
autant d'cquations que d'inconnues. II peut valoir la peine d'examiner ces so- 
lutions ad hoc. Pour autant que j'en puisse juger, elles combinent les seconds 
membres avec des coefficients qui sont toujours 1, ou —1. L'avantage est 
d'eviter les multiplications. L'inconvenient est de faire des choix ad hoc au vu 
des equations. Maintenant qu'il est d 'usage de faire des choix alcatoires, on 
peut se demander si des choix alcatoires de coefficients 1, et —1, suivant une 
loi ne dependant que du nombre d'equations et du nombre d'inconnues, ou 
plus simplement des choix aleatoires de 1 et — 1 avec la probabilite naturelle, 
ne donneraient pas des solutions presque aussi bonnes avec, par repetition du 
procede, une estimation assez correctc dc I'crrcur possible. Mcme s'il s'agit 
du centre de gravite, le choix alcatoire des coefficients 1 et — 1 a somme nuUe 
donne une estimation rapide de I'erreur possible, au lieu du calcul de la va- 
riance. 

Ce qui manquait a I'expose de Legendre etait justement I'estimation des 

erreurs possibles sur les quantites calculces. La liaison avec les probabilites 
allait etre faite par Gauss, et reprise par Laplace des 1810 dans son « memoire 
sur les integrales definies et leur application aux probabilites » ; les applications 
comprennent la theorie des erreurs, celle de la combinaison des observations 
(dont je viens de parler), et la methode des moindres carres, dont Laplace 
declare : 

« Cette methode, qui jusqu'a present ne presentait que l'avantage de four- 
nir sans aucun tatonnement les equations finales necessaires pour corriger les 
elements, donne en meme temps les corrections les plus precises. » 
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Bien avant Gauss, Laplace avait en main tons les elements pour lier les 
probabilites et la methode des moindres carres. L 'apparition de la fonction de 
Gauss et le calcul de son integrale sur la droite apparait chcz lui des 1774, 
dans son « mcmoirc sur la probabilitc des causes par les cvcncmcnts ». En 
effet, la loi d'une grandeur x quand on connait la loi des erreurs possibles 
et les resultats d'observation x\x",x"' est bien un calcul de probabilites des 
causes quand on connait les evenements. Laplace en donne la forme generale 

(p{x' — x)ip{x" — x)ip{x"' — x) 

quand </?(•) est la loi des erreurs. Le cas particulier ou if est la fonction de Gauss 
(Gauss est ne en 1777) est aussi un cas limite, et la methode des moindres 
carres avec expression des corrections aurait pu apparaitre des ce moment, 
mais ce n'a pas ete le cas. Laplace attribuait une grandc valcur a la thcorie des 
probabilites, il avait contribuc au calcul dc grandeurs a partir d'obscrvations 
en combinant des equations, mais le lien entre les deux lui a echappe. II a ete 
decouvert par Gauss avant que Laplace I'expose en 1810 puis dans son grand 
traite « Theorie analytique des probabilites ». 

La methode des moindres carres avec estimation des corrections etait-elle 
pratiquee avant d'etre parfaitement formalisee ? On pent se poser la question 
en lisant I'expose tres pedagogique qu'en fait Fourier dans deux memoires 
tardifs, de 1826 et 1828, dans le second desquels il mentionne en passant I'avoir 
appliquee lorsqu'il etait en Egypte, done avant 1800, pour la determination 
de la hauteur de la pyramide de Keops. Ces memoires ne contiennent pas de 
demonstrations, mais des enonces et des applications, et ils sont destines a des 
statisticiens. Depuis 1815, apres avoir ete demis de ses fonctions de prefet de 
I'lsere par Napoleon au retour de Pile d'Elbe, il s'etait installe a Paris et avait 
ete charge, grace au comte de Chabrol son ancien eleve, de la direction des 
services dc statistique du dcpartement de la Seine. 

Puis, apres son election tardive a I'Academie des sciences en 1817, il fai- 
sait figure de connaisseur en statistique. Son infiuence pent se mesurer par 
I'audience du prix Montyon de statistique dont il avait ete le rapporteur pour 
I'adoption commc prix de I'academie, et plus encore par le temoignage de 
Quetelet qui lui doit, dit-il, ses connaissances en probabilites. 

C'est le moment de parler de la transformation de Fourier. A la suite des 
series de Fourier, elle apparait clairement dans I'ouvrage principal de Fourier, 
la Theorie analytique de la chaleur, qui sous sa forme achevee a ete publiee 
en 1822. Elle etait dcja la dans le memoire de Fourier de 1811 qui lui a valu 
le Prix de I'Academie des sciences (ou plutot, les academies royales ayant ete 
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supprimees par la Revolution, de la premiere classe de I'lnstitut de Prance) 
et qui n'avait pas ete publie a I'epoque. On y trouve la formule d'inversion 

ct de nombreuses applications. Mais la encore, il y a dcs antecedents. Le plus 
important est la presentation par Clairaut de la transformation de Fourier dis- 
crete, celle dont il est question au debut de mon expose, dans son memoire de 
1754 sur I'orbite apparente du soleil ; naturellement, ni Clairaut ni personne ne 
pouvait soupgonner alors I'importance de la notion, et I'infiuence de Clairaut 
sur Fourier et ses successeurs semble avoir ete nuUe. 

Aujourd'hui le nom de Fourier est utilise principalement pour ce qui faisait 
le plus question a son epoque, les series de Fourier, les integrates de Fourier, la 
transformation de Fourier, et plus generalement I'analyse de Fourier. C'est a 
Rieniann que I'on doit I'appellation de series de Fourier pour les scries trigono- 
mctriques dont les coefficients sont donnes par les formules intcgrales connues 
comme formules de Fourier. Dans sa these sur les fonctions representables 
comme sommes de series trigonometriques Riemann fait une belle etude histo- 
rique du sujet, il en montre les antecedents, puis met en evidence le couplage 
entre les formules integrales donnant les coefficients et I'expression de la fonc- 
tion comme somme de serie trigonometrique, en declarant que Fourier est le 
premier a avoir compris la portee de ce couplage. Aujourd'hui nous pouvons 
exprimer ce couplage comme celui de I'analyse et de la synthese. Fourier avait 
insiste sur sa valeur universelle, et la un temperament est necessaire : la theo- 
rie ne pent s'appliquer qu'a des fonctions integrables, et alors la convergence 
de la serie pose probleme. 

Ainsi on doit considerer les formules de Fourier, celles qui expriment les 
coefficients a partir de la fonction et celle qui exprime la fonction a partir des 
coefficients, comme un programme plutot qu'un enonce. Comme enonce, il a 
amene a la recherche de contre-exemples instructifs. Mais comme programme, 
il est a I'origine de la fornialisation de la notion d'intcgrale, avec I'integrale de 
Riemann, puis celle de Lebesgue, puis celle de Denjoy, puis des distributions 
de Laurent Schwartz. II est aussi a I'origine des differents moyens d'exprimer 
la representation de la fonction, ou du signal, comme somme d'une serie tri- 
gonometrique : les procedcs de sommation au lieu de la seule consideration 
des sommes partielles, et les espaces fonctionnels dans lesquels les fonctions 
apparaissent comme des points. 

Parmi ces espaces fonctionnels le premier et le plus important est I'espace 
de Hilbert. On peut a volonte le considerer comme un f.'^, espace des suites de 
carre sommable, ou comme im L^, espace des fonctions de carrc intcgrable. 
On doit preciser quels sont les indices pour les suites, et de quelles fonctions 
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il s'agit. Le cas classique, qui fait I'objet du theoreme de Riesz-Fischer, est 
relatif a P{Z) et ^^(T), T etant le cercle R/Z. L'isomorphisme est assure par 
la transformation de Fourier : c'cst la correspondancc entre une fonction perio- 
diquc (dcfinic sur T) ct la suite dc scs coefficients de Fourier. Void des points 
de repere chronologiques. L'integrale de Lebesgue date de 1901. L'egalite de 
Parseval, ainsi nommee par Fatou, se trouve dans sa these de 1906 : elle ex- 
prime que la transformation de Fourier realise une isometrie de vers 
mais non que cette isometrie est surjective. En 1907, Frederic Riesz et Ernst 
Fischer, independamment ou plutot concurremment, out public deux series de 
notes aux Comptes rendus de I'Academie des sciences etablissant le theoreme 
et jetant les bases de ce que Fischer a appele une geometric des fonctions. La 
priorite n'est pas claire, et on pent aussi bien parler du theoreme de Fischer- 
Riesz. La clc pour passer de Parseval a Riesz-Fischer est le fait que est 
un espace complet. Et ce dernier fait montre la superiorite de l'integrale de 
Lebesgue sur celle de Riemann. 

Cependant raconter I'histoire ainsi est un renversement de son cours reel. 
En 1907, on ne disposait pas des symboles ni La notion d'espace norme 
n'existait pas encore, et celle d'espace norme complet encore moins. II fallait de 
longues periphrases pour parler de suites de Cauchy, et meme pour enoncer le 
theoreme. C'est un cas typique ou I'outil de la demonstration, le fait que est 
complet, est devenu le theoreme principal, et que pour I'enoncer en trois mots 
on a eu besoin de definir les termes aujourd'hui classiques de et d'espace 
complet : I'enonce sous cette forme n'apparait qu'en 1930 dans la Theorie des 
operations lineaires de Banach. En bref, le lemme est devenu theoreme et sa 
substance est passee dans les definitions. 

Le cadre general est apparu au cours des annees 1930. C'est celui des groupes 
abeliens localement compacts et de leur dualite. Ainsi L^(T) et ^^(Z) sont 
en dualite, de meme que L^(M) avec lui-mcme ou (^{'L/N'L) avec lui meme. 
La transformee de Fourier realise un isomorphisme isometrique entre L^{G) 
et L?{G), G et G etant deux groupes abeliens localement compacts en dualite. 
La notation i au heu de L est utilisee lorsque le groupe est discret. 

L'aspect geometrique est traduit aussi par la terminologie des bases or- 
thonormales et des systemes orthogonaux. Avec leur ancetre, la methode des 
moindres carres, et tous les developpements d'analyse lies aux systemes or- 
thogonaux, les suites de carre sommable out une place de choix dans I'univers 
mathematique. 
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Et cependant il y a plus naturel encore que et i'^. C'est et i^. En 
termes classiques, les transformees de Fourier fournissent les series de Fourier- 
Lebesgue (les coefficients sont donnes par les formules de Fourier avec comme 
integrale I'integrale de Lebesgue), appelees aussi simplement series de Fourier, 
et d 'autre part les fonctions continues sommes de series trigonometriques ab- 
solument convergentes, qui constituent I'algebre de Wiener. Le point essentiel, 
mis en evidence par Wiener, est que ces espaces et sont des algebres pour 
la convolution. Par transformation de Fourier, convolution et multiplication 
s'echangent. La convolution se trouve partout en analyse, mais en tant que 
notion c'est I'une des dernieres nees : le terme de convolution ne s'est impose 
que dans la seconde moitie du 20'' siecle. Au moins Wiener en avait fait sous le 
nom de Faltung un merveilleux usage, qui a servi de modele a plusieurs egards 
a la theorie des distributions de Schwartz. 

L'algebre de Wiener a ete un grand sujet de I'analyse au cours des annees 
1950 et 1960, et j'ai collabore a cette epoque a la decouverte de ses secrets. 
Puis d'autres sujets sont passes au premier plan. II est remarquable pour moi 
que les problemes actuels de theorie du signal et de traitement des donnees 
donnent un nouveau depart a ce secteur de I'analyse de Fourier. 

Wiener est connu a bien des titres. Le pense ici au mouvement brownien 
et aux series de Fourier a coefficients aleatoires, un autre sujet sur lequel j'ai 
travaille dans la tradition de Paley, Wiener et Zygmund. Le theoreme de Can- 
des, Romberg et Tao par lequel j'ai debute cet expose marie la transformation 
de Fourier discrete, I'espace £^ et les probabilites. C'est done de plusieurs ma- 
nieres un heritage de Wiener. 

II est bon de conclure en revenant a ce theoreme pour en indiquer le contexte 
et d'eventuels prolongements. 

L'article |CRT06] donne en introduction I'expose d'une experience nume- 
rique qui donne d'excellents resultats et qui viole I'usage etabli par le theoreme 
d'echantillonnage de Shannon. II s'agit d'une figure appelee le « fantome de 
Logan-Shepp », constituee par quelques plages en forme d 'ellipses qui em- 
pietent les unes sur les autres, et qui est un test classique pour les methodes 
d 'analyse du signal. On regarde sa transformee de Fourier sur 22 demi-droites 
issues de I'origine et regulierement reparties, puis les valeurs qu'elle prend sur 
une vingtaine de points sur chaque demi-droite, au total 512 points. 

On reconstitue exactement la figure par un procede d'analyse convexe a 
partir des valeurs prises par la transformee de Fourier sur ces points, qui sont 
espaces les uns des autres d'au moins 7r/22, alors que le theoreme d'echan- 
tillonnage exigerait pour etre applique que I'espacement des points voisins 
soit 27: /N avec N = 1000 environ. Cet exemple n'est pas encore parfaitement 
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compris, mais il repose sur la recherche du minimum dans I'algebre de Wiener 
des normes des fonctions qui extrapolent la restriction a un ensemble conve- 
nable de la transformee de Fourier du signal considere. Le cadre des groupes 
cycliques d'ordre A'' pour G et G se prete bien a I'etude theorique, et I'echan- 
tillonnage aleatoire aboutit a I'enonce que j'ai indique. 

L'echantillonnage aleatoire est essentiel dans la demonstration, et il est 
egalement tres interessant en pratique. En experimentant son utilisation on 
constate dans des cas simples que le choix de I'effectif de I'echantillon, pent 
etre rendu notablement plus petit que la valeur fixee par le theoreme. L 'article 
signale les resultats d'experiences allant dans ce sens. 

Quel est le champ d 'application de ces methodes ? Au contraire de la me- 
thode des moindres carres, qui consistait a gerer un ensemble de donnees 
surabondantes, il s'agit ici de restreindre les observations au minimum indis- 
pensable, un peu plus de T si Ton salt que le signal est constitue de T pics au 
plus. 

L'utilisation de I'analyse convexe, qui rejoint ce que Fourier appelait I'ana- 
lyse indeterminee, c'est a dire I'etude des inegalites et leur interpretation geo- 
metrique, s'est introduit progressivement dans ce domaine. Les normes £^ se 
sont averees d'usage commode, lorsque la nature des problemes suggerait la 
norme i^, c'est-a-dire le nombre de termes non nuls. Un exemple typique se 
trouve dans I'article de Donoho et Huo de 2001 dont la motivation est la re- 
construction a partir d'un developpement mixte en ondelettes et en Fourier, 
ou en ondelettes et en bandelettes, d'une image comprenant des plages et des 
contours. Le cadre general est celui deja dresse en 1993 par Mallat et Zhang : 
on dispose d'un « dictionnaire » constitue de plusieurs bases orthonormales, 
et il s'agit de I'utiliser pour obtenir de meilleurs resultats qu'avec une seule 
base. Dans Particle [DHOl] , I'etude theorique de depart concerne les groupes 
cycliques d'ordre A'", G et G. Toute fonction definie sur G pent se traduire 
par Fourier en une fonction definie sur G, et reciproquement. II y a done 
plusieurs representations possibles pour la meme fonction, en exprimant une 
partie sur G et une partie sur G. Une representation economique utilise le 
minimum de points de G et de G, soit T points de G et 17 points de G ; sa 
norme dans est alors T -|- J7 ; et on pent esperer obtenir la meilleure repre- 
sentation en minimisant la norme i^. C'est bien le cas lorsque T + ^ < ^/N, 
ce ne I'est plus lorsqu'il y a egalite (description de I'exemple ? c'est la fonc- 
tion indicatrice du sous-groupe d'ordre n quand N = ■n?, qui, a un facteur 
multiplicatif pres, est sa propre transformee de Fourier : le minimum de la 
norme est atteint de deux manieres). Mais chercher le minimum dans 
est impraticable. Le minimum dans est praticable et donne le bon resultat 
quand T-FO < \^ . 
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PARTIE II 

RESTRICTIONS ET PROLONGEMENTS 
DANS L'ALGEBRE DE WIENER 



Introduction 

Les travaux sur rechantillonnage parcimonieux (compressed sensing) en 
theorie du signal ont mis en relief un tres interessant probleme d'analyse de 
Fourier. En me referant aux notations d'Emmanuel Candes dans |Can06] . je 
le presenterai de la maniere suivante. 

G est un groupe abelien localement compact et G son groupe dual. Soit (5) 
une partie de L^(G), et K une partie de G. A toute x E L^{G) on associe 
sa transformee de Fourier x, qui appartient a I'algebre de Wiener A(G) (= 
FL^{G))^ et la restriction x a x\k- Les y G L^{G) tels que y\K = x\k 
constituent un convexe ferme dans L^{G), et les y, prolongements de x\k 
dans A(G), constituent un convexe ferme dans A{G). II s'agit de donner une 
condition sur (S) et K garantissant que, si x G (S), x est le prolongement 
minimal unique de x\k dans A{G) : 

(1) ll^oll^(g) = inf{||y||^(g) I y\K = x\k} xo = X. 

Dans |Can06j et dans les travaux qui precedent, de Candes, Romberg et Tao 
|CRT06| et anterieurement, sur un sujet voisin, de Donoho et Huo [DHOl] , 
G et G sont le groupe cyclique TLjNTL^ N etant un entier tres grand. Le groupe 
G est une representation discrete des temps et le groupe G des frequences. En 
designant par S le support du signal x, (S) est traduit par une condition sur 
S, en I'occurrence dans [Can06j et |CRT06] que le cardinal de S, \S\, est 



majore par un nombre T petit par rapport a : 

(2) \S\ = \\x\\io T <^ N 

(la notation || \\io a justement (2) pour definition). K est un ensemble de 
frequences, de cardinal \K\. Le theoreme principal de Can06 et |CRT06] 



est que, si I'on fixe \K\ et que I'on choisit K aleatoirement avec la probabilite 
naturelle, on retrouve le signal x par le procede (1) avec une probabilite voisine 
de 1 si \K\ est assez grand par rapport a Tlog A^, a savoir que 

(3) I^I^CriogiV 
avec C ^ 23(M + 1) entraine 

(4) P{{1)) ^ 1 - cN-^ , 
ou c ne depend que du signal. 
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Ainsi, si Ton sait que le signal x est porte par T points, avec T ^ A^, on 
est loin d 'avoir besoin des N coefficients de Fourier de x pour le reconsti- 
tuer : en prenant C = 100 dans (3), \K\ coefficients choisis au hasard suffisent 
pratiquement. 

Les titre de |CRT06] et |DH01] mettent en avant ce que les auteurs ap- 
pellent « uncertainty principles », c'est-a-dire le fait (incontestablement cer- 
tain) que les supports de x et de x ne peuvent pas etre petits simultanement. 
Les relations de Heisenberg, d'ou est tire le vocable de « principe d'incerti- 
tude », expriment simplement que Ton ne pent pas localiser simultanement 
un signal en temps et en frequence. En physique, elles signifient que la re- 
presentation des particules par des points materiels est inadequate a I'echelle 
quantique. En mathematiques, elles constituent un theoreme de localisation 
sur transformees de Fourier. Pour G = Z/NZ, G est une autre copie de Z/NZ, 
et le theoreme de localisation s'exprime ainsi : 

(5) \\x\\e° ■ \\^\\e° ^ 

La demonstration, que j'indique en appendice, donne des informations supple- 
mentaires en fonction des proprietes arithmetiques de A^. On exprime souvent 
le theoreme de localisation dans les Z/AZ par une consequence de (5), a savoir 

(6) \\x\\io + \\x\\io ^ 2\/A. 

L'egalite dans (6) a lieu lorsque A est un carre parfait et que les supports de x 
et de x sont des translates du sous-groupe cyclique d'ordre \/A du groupe 
Z/AZ, en particulier quand x et x sont le « peigne de Dirac » d'ordre \/A, 
divise par \/A- 

Dans [DHOl , Donoho et Huo tirent parti de cet exemple comme limitation 
naturelle pour la question qu'ils posent de la representation d'un signal la plus 
economique possible en temps et en frequences, c'est-a-dire I'expression de x, 
fonction definie sur G = Z/AZ, sous la forme 

(7) X = y + z, \\y\\e° + 112^11^° minimum. 
Quand A est un carre parfait et que 

(8) W^Wi" = = V^, 

on pent choisir x = y ou x = z dans (7) : la representation la plus economique 
n'est pas definie. Donoho et Huo montrent qu'elle est bien definie lorsque x 
est de la forme 



(9) 
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et que, sous I'hypothese additionnelle 

(10) x = y + z, llyll^o + pll^o < ^\/iV, 
on obtient la solution (7) en remplaQant la condition 

(11) l|y||^° + 11^11^° minimum, 
pratiquement ingerable, par la condition 

(12) \\y\\e^ + ll^ll^i minimum, 

qui se prete au calcul. 

Candes, Romberg et Tao utilisent leur theoreme sur le prolongement mi- 
nimal (formules (3) et (4)) pour proposer une amelioration considerable du 
theoreme de localisation (6) sous la forme d'une proposition hautement pro- 
bable : si I'on choisit au hasard une partie S de G = 'LjN'L et une partie S' 
de G telles que 

(13) |S| + |5'|<7 ^ 



alors, avec une probabilite aussi voisine de 1 que Ton veut par choix de 7, il 
n'existe aucune fonction x definie sur G a support dans S telle que x ait son 
support dans S' . 

Par cet aperQu tres partiel des travaux de Donoho et Huo de 2001 et de 
Candes, Romberg et Tao de 2006, on voit la portee en theorie du signal de 
la methode de minimisation dans ^. Cette methode n'est pas nouvelle, et 
son histoire est retracee brievement dans |CRT06] p. 494. Les initiateurs, au 
cours des annees 1980, semblent avoir ete des geophysiciens. La minimisation 
des normes £^ est une question d'analyse convexe dont on pent trouver les 
prolegomenes dans ce que Fourier appelait I'analyse indeterminee |Fou90] . La 
reduction de ce sujet d'analyse convexe a des techniques de programmation 
lineaire a ete exposee en 1998 par Chen, Donoho et Saunders dans |CDS98 



La minimisation de normes ^ est aujourd'hui un outil important de la theorie 
du signal. 

Je vais m'attacher aux conditions de validite de (1) dans un contexte tres 
different : la reconstruction d'une fonction a spectre lacunaire a partir de sa 
restriction a un intervalle par prolongement minimal dans I'algebre de Wiener. 
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Reconstruction des fonctions a spectre lacunaire 

Dans cette section on s'ecarte de I'analyse de Fourier discrete. Le groupe G 
sera Z'^ (z^ = 1, 2, . . . ) et G sera T*^. On s'interesse aux fonctions x G ^(T'^), 

(14) x{t) = ^ x{n)e{n ■ t), \\x\\a = ^ |2;(n)| < oo 

{t = {ti, . . . ,ty), n = (ni, . . . ,nj,), n ■ t = riiti + • • • n^ty, e{u) = exp(27rm)) 
dont le spectre S est lacunaire au sens que deux points distincts de S sont 
toujours a une distance ^ d : 

(15) inf Is' - s| ^ d 

s'£S 

I • I represente la distance euclidienne, et le premier membre de (15) s'appelle 
le « pas » de S. Ici K sera une boule fermee de rayon r dans T*^ (0 < r < 1/2). 
Le probleme pose dans I'introduction est de donner une condition sur d et r 
(z^ etant donne) garantissant que x est le prolongement minimal unique dans 
A(T'^) de sa restriction a K, x\k, c'est-a-dire que la formule (1) est valable. 

Pour faciliter la comprehension, limitons-nous d'abord au cas z/ = 1. Alors K 
est un intervalle, et comme la solution du probleme est invariante par les 
translations sur T, on pent choisir K = [— r, r]. Ecrivons (14) sous la forme 

(16) x=^3;„e„, \\x\\A = Y\xn\ < oo- 

n€S n 

(1) signifie que, quelle que soit z G A{T) nulle sur K et ^ 0, soit 

(17) z=^2;„e„, 0<^|2;„|<oo, z\k = ^, 

raeZ n 

on a 

(18) Y\Xn + Zn\ + Y\Zn\>Y\Xn\■ 
n&S nfS n€S 

Dire que (18) a lieu pour toute suite (x„) G ^^(S), c'est dire que 

(on a utilise sup^j^dx^l — \xn + Zn\) = \zn\ ; c'est le meme calcul qu'en [CRT06I 
p. 494). On cherche done une condition sur d et r pour que (17) implique (19). 

II y a bien des fagons d'aborder le probleme. Voici celle qui m'a donne les 
meilleurs resultats (loin sans doute d'etre optimaux). On multiplie z, qui est 
nulle sur K, par une fonction (p qui est presque entierement portee par K. Le 
resultat est voisin de O, et il en est de meme de sa transformee de Fourier, 
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convolution de z et de Si est presque entierement portee par I'intervalle 
[— d, d], la quasi-nuUite de la convolution 2; * est capable d'entramer (19). 
II est commode de partir de la fonction de Gauss 

+2 



7(0 



(20) 
de poser 

(21) 7a(t) = a^{a£) 
et de reduire 7a modulo 1 sous la forme 

(22) 7;(«) = > -in\t-m) 



{t e M), 

(a » 1/r) 



^7(ra/a) e„(i). 



7:(i) = E 

C'est 7* qui nous servira de fonction ip. Comme elle est paire, et decroissante 
et convexe sur I'intervalle [r, sa restriction a T \ if a une norme ^ 27* (r) 
dans I'algebre quotient ^(T \ K), d'ou resulte que 

(23) ll^7alU(T) ^ 2||^|U(T)7a('')■ 

Pour chaque s E S, integrons sur T le produit £7*6-5 : 



(24) 



E^"7 



,n 



= Zs + Y^ Znl 



,n 



Multiplions par Eg = \zs\/zs et faisons la somme sur S : 



(25) 



Zlle-s = E l^^l + E ^sZnl{ 
sES n&'Z,s&S,n^s 



n 



'-). 



Dans cette derniere somme, distinguons les n € 5*, que nous appellerons s', et 
les n ^ S, et remarquons que le premier membre de (25) est majore en valeur 
absolue par celui de (23). On obtient 



(26) 



avec 



(27) 



E \zs\ + E Ms')zs' + E ^(")^n ^ 2||£|U(T)7:(r-) 
ses s'es nis 



E ^si{ 

sG5\{s'} 



ses 



et (26) entraine 



(28) E \^n\mn)\ + 27:(r)) ^ E 1^^1(1 " 1^(^)1 " 27:(r)). 
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En vue d'obtenir (19), supposons 

1 



(29) 



< a < d et a ^ 10 



et majorons |i?(n)|, 7Q(r) et |yl(s)|. 

Quand n ^ S, soit h la distance de n a S". Alors 



(30) 

Pour s G iS, 
(31) 

Enfin 
(32) 

Supposons 

(33) 

alors 

(34) 

Supposons 

(35) 

alors 

(36) 



et (19) est bien etablie. 
L'une des conditions 

(37) 

ou 
(38) 



\B{n)\ ^ ^{h/a) + 7((c/ - h)/a) + 2 J2 lijd/a) 

3=1 

oo 

^j{l/a)+j{d/2a)+2j2l{jd/a) 
^ 1 - (37r/4a2) + 2exp ( - nd'^ /Aa^) . 



\A{s)\ ^ 2^7(id/a) ^ 2exp ( - Trd^Aa^). 



27* (r) ^ 3aexp(-7ra^r^). 



r > - \/log a ■ 
a 



27: (0 ^ 



TT 



d > 2a-\/loga : 



TT 



1^(^)1 < ^ 

|S(n)|^l- 



27r 
4^ 



^ 5^ 1 
a > - log - 
r r 



^^<10^ 



r > - log d (d > 10) 
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entrame done, en choisissant a = y^d/r, ce qui implique (29), que pour toutes 
les fonctions z verifiant (17) on a (19). 

L 'etude faite en dimension 1 se traduit sans peine en dimension v > 1, au 
prix d'unc augmentation dcs constantcs. 

On prend pour K la boulc |t| < r ct on rcmplace 7, 7^ et 7* par 

(39) 7^0 = e-"l*l' {t G R^) 

(40) 7„,,(t) = a^jiat) (a > ^) 

7a,i.(*) = 7a,:/ (* " "^) 

(41) 

= 7i/("/a)e(ra-i) 

La formule (28) reste valable en remplagant 27* (r) par ||7a ^11^(1" 
et A{s) et B{n) par Ai,{s) et B^{n), definis comme en (27) avec ji, au lieu 
de 7. La majoration de ||7a,;y|U(T''\_ft') peut se faire ici en estimant les derivees 
d'ordre ^ de 7* sur \ K, et au lieu de (32) on a 

(42) ||7a,JU(T-^x/f) ^ C'z/a''exp(-7raV), < 

La majoration de |Biy(n)| se fait comme en (30), en distinguant le point s le 

plus voisin de n (ou I'un dc ccs points s'il y en a plusicurs) et les points de S 
dans la boulc dc centre s ct dc rayon |, dont le nombres est ^ 2^. Celle de 
|Aiy(s)| est analogue, en plus simple. On obtient 



(43) 



\Bin)\ ^l-^ + 2^+iexp(-.^) 

|yl(.)| ^2-+iexp(-vr^), 



d'ou finalement I'enonce suivant. 

Theoreme. — Soit G = T'^{u = 1,2,...). Pour que toute fonction x G A{G), 
dont le pas du spectre est d (d'^ 10) s'obtienne d partir de sa restriction x\k 
d une boulc K de rayon r comme prolongement minimal de x\k dans A(G), 
il suffit que I'une des conditions suivantes soit realisee : 

5 

(44) ^ ^ d ^og{i'd) 



(45) d>^ ^log(l/r) + 1 
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(Pour = 1, ce sont les formules (38 et (37)). 

La constante 5 peut etre abaissee quand d est grand, et doit etre ajustee 
pour couvrir tous les cas. 

Discussion et commentaires 

Je commencerai par le cas 1^ = 1. 

1. On peut transcrire le theoreme dans le contexte donne dans I'introduc- 
tion, ou G et G sont le groupe cyclique Z/iVZ, N tres grand, G etant une 
representation discrete des temps ct G dcs frequences. Alors x est un signal 
porte par des temps scparcs par dcs intcrvallcs dc longueurs ^ d, ct K est 
une bande de frequences de largeur ^ 2rN (si Ton represente G par les ra- 
cines AT-iemes de I'unite, K est la partie de G contenue dans un arc du cercle 
unite) . On obtient I'enonce que voici : 

Si le support S du signal x (porte par Z/iVZ) a un pas ^ d, c'est-a-dire si 

(46) inf \s' -s\^d (d ^ 10) 

s'es 

on peut obtenir x a partir de frequences consecutives (dans Z/NZ) par le 
procede indique en (1), des que 

N , 

(47) \K\ ^W — 

Lc dcfaut dc cct enonce par rapport a celui de Candes, Romberg et Tao est 
dans riiypothcsc faitc sur S : (46) implique \ S\ ^ N/d, mais jS"] ^ T n'implique 
rien sur d. Cependant, si T est donne et si S s'obtient en disposant T points 
au hasard sur Z/iVZ, on a N/d > lO/T^ avec une probabilite > 0,95, done il 
est probable que 

(48) \K\ ^ lOOT^v/bgT 

suffit pour que (1) soit valable. L'avantage est qu'il s'agit de frequences conse- 
cutives, c'est-a-dire d'une bande de frequences. 
Cela serait a comparer aux resultats connus. 

2. On a beaucoup etudie les series trigonometriques lacunaires et les fonc- 
tions qu'elles representent : 

(49) m = J2 ^nenit) = Yl ^ne™ {t e T). 

neS neS 
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C'est de ces series et fonctions que nous nous sommes occupes a partir de (16), 
sous la condition X^l^nl < oo. On pent les etudier sous d'autres conditions, 
en particulier < cxd. 

En gros, si S est lacunaire, beaucoup de proprietes de x se lisent sur les 
restrictions x\k, ou K est un ensemble convenable. Je me limiterai ici au cas 
des intervalles. 

Premiere propriete : la nuUite. Quand est-il vrai, sous I'hypothese (49), que 
x\k = entraine x = ? Autrement dit, que x soit bien defini par les valeurs 
qu'elle prend sur K ? Cela se traduit lorsque K = [— r, r] par une question sur 
les fonctions entieres de type exponentiel ^ 7r(l — r) bornees (ou L^, ou L^...) 
sur la droite reelle, X(z) (z G C) : a quelle condition leur nullite sur Z \ 5" 
entraine-t-elle leur nullite ? Plus simplement, S etant donnee, quelle est la 
borne inferieure des r tels qu'il en soit ainsi ? La reponse est fournie par la 
theorie de Beurling et Malliavin [BM62|, IKoo98 . Associons a toute suite 
A C Z sa fonction de decompte n{t), dont la variation sur un intervalle egale 
de nombre de points de A sur cet intervalle, et convenons de dire que A admet 
A comme densite-BM si 

(50) / Mi^,, < ^. 



l + t2 

La densite-BM interieure de S se definit naturellement comme la borne su- 
perieure des densites-BM des suites A contenues dans S. C'est justement la 
borne inferieure dans \K\ (longueur de K) tels que x\k = entraine x = 
{K=[-r,r]). 

En particulier, pour que x soit bien definie par les valeurs qu'elle prend sur 
n'importe quel intervalle, si petit soit-il, il faut et il suffit que S ne contienne 
aucune suite A ayant une densite-BM > 0. II n'est pas necessaire que S soit 
tres dispersee ; il suffit qu'elle ait de grandes lacunes. 

Seconde propriete : I'appartenance a I'algebre de Wiener. Rappelons que 
A{T) = T£^{Z) et que A{K) est I'algebre quotient de A{T) par I'ideal constitue 
par les fonctions nuUes sur K, autrement dit, I'algebre des restrictions a K des 
fonctions € yl(T). On considere les fonctions x de (49) avec x G Quelle 
est la borne inferieure des longueurs des intervalles K tels que x\k G ^(-^) 
implique x £ A(T) ? Reponse : c'est la densite uniforme exterieure de S, c'est- 
a-dire la borne inferieure des densites D des suites A telles que 

(51) n{t) = Dt + 0{1) {t^±oo) 

et qui contiennent S ; on exprime (51) en disant que A a D pour densite 
uniforme [Kah70 . 
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Variante : I'appartenance aux espaces A^'^iT) = J^i^'^'^Z). La definition de 

^''■(Z) (1 a ^ oo, r G M) est 

{{x I EnGZ knr(l + <oo| sil^a<oo, 

X I SUp^g^ |2;n|(l + I"-!)*" < CXD ^ si a = CXD. 

En general, A°'''^(T) est un espace de distributions de Schwartz, dont on sait 
definir les restrictions a des intervalles ouverts J; ces restrictions constituent 
I'espace A°'''''{J). Question : quelle est la borne inferieure des r tels que x\j G 
A°''^(J) (J =] — r,r[) entraine x G A"''^(T) ? Reponse : elle est la meme que 
precedemment, independante de a et r. 

Pour r = et a = 1, on retrouve ^(T). Pour r = et a = 2, on trouve 
L'^(T), et c'est I'origine de la theorie. Sans se borner a des A C Z mais en 
considerant des A reelles « regulieres » au sens que leur pas 

inf IV- A| 

aga 
A'eA 

est > 0, Paley et Wiener ont appele pseudo-periode attachee a A la borne 
inferieure des ^ > tels que les normes dans L'^{J) des fonctions 

(53) /(t)=^XAe^^* 

soient equivalentes pour tons les intervalles J de longueur £. Au facteur 27r pres, 
la valeur de la pseudo-periode est la densite uniforme exterieure de A |Kah62 . 

Troisieme propriete : la continuite. C'est la question la plus naturelle, mais 
la plus difficile. Je renvoie a I'etude d'Yves Meyer |Mey73] pour les meilleurs 
resultats connus. 

Quatrieme propriete : I'analycite. Sous une forme un peu differente, c'est 
la premiere question qui a ete examinee (Fabry en 1896, Polya en 1942), mais 
la derniere ayant eu sa solution completement explicitee |Koo98| II p. 52] 
|KLR98] p. 199]. Quelle est la borne inferieure des longueurs |J| telles que, 
si X est analytique sur J, x est analytique partout sur T? C'est la densite 
maximale de S, a savoir la densite exterieure quand, au lieu de (50) et de (51), 
la densite D d'une suite A est definie par 

(54) n{t)=Dt + o{t) {t^±oo). 

Polya a montre que la densite maximale est atteinte : il existe une suite conte- 
nant S dont la densite (au sens de (54)) est la densite maximale de S. 



ANALYSE ET SYNTHESE HARMONIQUES 



25 



Lorsque le pas de S est ^ d, toutes les densites exterieures de S sont ^ 1 /d. 
On pourrait done esperer reeonstruire x a partir de x\k {K = [—r,r\) quand 
\K\ = 2r > 1/d. Nous avons vu (formule (44)) que c'est possible quand \K\ > 
^^/^ogd. On pent ameliorer le resultat en diminuant le second membre, mais 
pas au dessous de 1/d. 

Aucune densite n'est adaptee a la reconstruction d'une fonction x appar- 
tenant a A(T) a partir d'une restriction x\k par prolongement minimal dans 
A(T). En efFet, soit z £ A(T) une fonction nuUe sur K. Soit x une somme 
partielle de sa serie de Fourier telle que H^ — x||yi(Tr) < l|2|U(T) ! alors le pro- 
cede (1) echoue pour x, et le support S de x est fini, done de densite nuUe. 
L'hypothese naturelle pour cette reconstruction est bien que la pas de S soit 
convenablement minore. Le theoreme de Candes, Romberg et Tao a le grand 
merite de prendre pour hypothese une majoration de IS"! sans rien supposer de 
son pas ; mais il faut alors renoncer a prendre pour K un intervalle, et choisir 
pour K un ensemble fini aleatoire. 

Quittons le cas v = 1. 

3. Pour 1/ ^ 2, on est loin d'avoir des reponses aussi precises. Par exemple, 
on ne salt pas pour quelles valeurs de r l'hypothese 

(55) inf \s' — s| ^ d, 

s'&S 
s'^s 

I I representant la norme euclidienne, entraine I'equivalence des normes 
II^IIl2(_b) pour toutes les boules B de rayon r. La relation entre la geome- 
trie de S et ceUe des K associes est un probleme interessant, amorce dans 
|Kah62 . 

Appendice 1. Sur le principe de localisation de Tao 

G et G etant deux groupes cycliques d'ordre en dualite, x une fonction 
definie sur G a valeurs complexes et x sa transformee de Fourier, il s'agit 
d'etablir que 

(56) ■ WxWt" ^ N. 

Void la preuve de Tao, ingenieuse et facile. 

Soit S le support de x, et IS*! = \\x\\iio = s ^ 0. Rappelons que 

x{t) = x„e2™*/^. 
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Supposons x{ti) = x{t2) = ••• = x{ts) = avec ti,t2, ■ ■ ■ ,ts en progression 
arithmetique, tons distincts modulo N. Les x„ (n G S) sont solutions d'un sys- 
teme lineaire homogene dont le determinant, du type Van der Monde, est 7^ 0, 
done ils sont tons nuls, contrairement a I'hypothese. Ainsi, sous I'hypothese 
llxll^o = s > 0, X ne pent pas s'annuler en s points formant une progression 
arithmetique. 

Les intervalles d'entiers contigus au support S" de x ont done au plus s 
elements, ce qui entraine s\S\ ^ N, c'est-a-dire (56). 

On a egalite dans (56) lorsque et sont des translates de sous-groupes 
cy cliques d'ordre s et s, de G et G, avec ss = N. 



Appendice 2. Reconstruction d'un signal par la methode 
de Candes, Romberg et Tao 
(variantes deterministe et probabiliste) 

Le signal est une fonction x{t), t G G = Z^v = Z/NZ, portee par une 
partie S de G, de cardinal T. Sa transformee de Fourier est 

x{t) = ^Y1 x{t)e{-ujt/N), w G G = Zjv, e{u) = e^™. 

On veut reconstituer le signal en n'utilisant qu'une partie O des frequences lo, 
par le procede que voici : 

(i) X est le prolongement minimal de x\q dans A{G) (= J-L}{G)). 
Quand est-ce possible ? On va se donner des conditions suffisantes, emboi- 

tees. J'ecrirai indifFeremment x{t) ou xt- (0) signifie que pour toute z nulle 
sur O et ^ 0, 

^\xt + zt\+ ^ \zt\>^\xt\ 
teS tGG\5 tes 

et il suffit pour cela que 

Done il suffit pour avoir (0) que 

(ii) Pour toute z nulle sur et ^ 0, J2tGG^S \^t\ > Z^tes \ ^t\- 

Pour traduire (jn]), introduisons At {t £ G) telle que |At| = 1 et XtZt = \zt\. 
Soit P{t) (t £ G) tel que P soit porte par il. Alors 

Pz= fPz = 0, 
G Jg 
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soit 
(1') 

Supposons que 

(20 

Alors 



+ 



Gx5 



\P{t) - Xt\ < 1/2 
\Pit)\ < 1/2 




Pz = 0. 

quand t £ S, 
quand t £ G \ S. 

d'ou 



/ Pz 




"■lis 


Js 





et, compte-tenu de on a 



\z\ > 



G\5 



c'est-a-dire la conclusion de (jii]). II suffit done, pour avoir (juD, que 

(iii) quels que soient les Xt {t £ G) de module 1, il existe P{t) {t £ G) tel 

que P soit porte par O et qu'on ait (f2^ . 

Reste a construire P. Pour cela, introduisons le noyau 

(3') Kit) = e{ujt/N) {t £ G) 

et posons 

K{t - t') 



P{t) 



t'£S 

1 



K(0) 



^ K(0) ^ ^' 

Lu£Q t'es 



' e 



-Lot'/N) e{ujt/N). 



Ainsi, P est porte par J7. Supposons maintenant que 
(iv) le noyau K{t) defini par ([s]) verifie 

\K(t)\ < ^ K{0) quand t / (rappelons que T = \S\). 

Alors (|2^ est verifiee done dm]) a lieu, done (|ii|), done (jl]). La condition (jiv|) 
entraine (jl]) pour tous les signaux x{t) portes par T points. 

Ce resultat faeile est la variante deterministe du theoreme CRT. Je eompa- 
rerai plus tard la portee du theoreme et des variantes. 

La reeherehe des noyaux (p^ ayant la propriete (jiv|) est interessante en elle- 
meme. Les noyaux s'appellent idempotents paree qu'ils eoineident avee leur 
earre de eonvolution ; la propriete (|Iv]l traduit une eertaine parente entre O 
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et le groupe G entier. Pour aborder cette recherche, nous allons emprunter 
I'approche de Candes, Romberg et Tao. 

Selon |CRT06 , choisissons maintenant pour Q la partie de G definie par 
une selection aleatoire de parametre r : pour chaque n G G, I'evenement 
{n G 17} a pour probabilite r, et ces evenements sont independants. En posant 

les Xn sont des variables de Bernoulli de parametre r, independantes : 

distribution de = -6(1, r). 
Le cardinal de 0, est J2-^n, qui est voisin de Nt : 

distribution de \^l\ = B{N,t), 

et quand est fixe, la distribution de I'ensemble aleatoire Q est uniforme 
sur I'ensemble des parties de G a |0| elements. Dans la suite, ce qui sera dit 
de I'ensemble aleatoire Q sera egalement valable pour un ensemble de [Nt] 
elements distribues au hasard sur G selon la probabilite naturelle. 

Precisons ce point. Pratiquement, et Nt sont grands et on pent se repre- 
senter B(N, r) comme une distribution normale M{Nt, Nt{1 — r)) (esperance 
Nt, variance Nt{1 — r)). Si, comme dans la suite, tN = G(T) logN, on aura 

(4') -riV| > AVlogiV) = 0(iV-5^'/^(^)) (iV ^ oo). 

Comme exercice, cherchons a demontrer cette formule. Soit u > ; on a 

(5') S(e"^") = 1 - r + re" < exp(r(e" - 1)). 

Si de plus n < 1, 

P{ Ef Xn-TN- AVl^ > 0) < S(exp (n(Ef X^-tN- A^b^))) 

< exp(iVrii^ - AulogiV). 
Le minimum du second membre est atteint quand u = {XlogN)/2NT et vaut 

A2 log^ATs 



exp 



/ _ A- log iV \ 
V 4/Vt / ■ 



ANt 

En remplagant J2i Nn — tN par Nn + tN et en ajoutant les probabilites 

obtenues, on a 

A2 log^AT, 



^(1 El Nn - tN\ > X^/hiN) < 2 exp 



4Nt 
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c'est-a-dire, quand tN = C(T) log A^, 

(6') P{\M - tN\ > AVl^) = 2iV-(^'i°s'^)/4^\ 



L'exposant 1/4 dans ([64) au lieu de 1/2 dans (j4^ tient a la majoration trop 
genereuse de — 1 par u + . 

Nous identifierons desormais G* a Z^r. Au lieu de ([s]), nous avons 

= Xne{nt/N) 

et en particulier i^(0) = X^^n- Etant donne N et T, on va chercher a deter- 
miner r de fagon que ()iv|) soit tres probable. 
Le contraire de ()iv|l est 

3tGG\{0}, |K(t)|^±i^(0). 

Choisissons un entier z/ ^ 3. Pour chaque z G C \ {0} il existe un 93 E 
{2fc7r/zv I A; = 1, 2, . . . , z/} tel que Re ze"**^ ^ |z| cos vr/i/. Posons a = cos vr/z^. 
Le contraire de ()iv|) implique 

Fixons t et posons 

y = ReK(t)e-^-^K(0) 

= ^ X„(cos(2^nt/iV -'/')- ^ 

et cherchons a majorer P{Y ^ 0). On a 

P(y ^ 0) < E'exp-uy pour tout li > 0, 

et, utilisant (Is^ . 

£'expuy= JJ^ exp |^m( cos(27rnt/A^ — 99) — — 

= exp exp cos(27rnt/A^ — (^) — — )^ — 1 



exp 



00 t. 



^ E E ^ ( cos(27rnt/iV " " ^ 



E ^ E ( cos(27rnt/7V - 9^) " 

fc=l ' n£Z]\i 
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Supposons impair. Comme t ^ 0, done aussi 2i 7^ 0, la derniere somme vaut 

.uand.^1, 



h N 



quand k = 2 



et elle est toujours majoree par A^2'^, done 



EexpuY ^ exprA^ 



a /I ^ 
2r^^ V2 ^ 4r2. 



^ fc=3 



Pour minimiser la somme des deux premiers termes du crochet, prenons 



u 



a /I a 
2r V2 ~^ 4r2 



On obtient 

E exp uY ^ exp tN 



a2 /I 



+ 



^ exp 



rA^ 



8r2V2 4r2 

2 \ 



a2 ^-1 ^ ^^/2axfc 



) +1: 



fc=3 



4r2 + 2 r3 / 



done 



P(non(jivD) < A^z^exp 



tN 



a" 2 \ 

4r2 + 2 ~ T^) 



Si Ton choisit 

rA^ = 4C(r2 + 1) log Af, C = l + 5>1, 

puis 1/ assez grand, done a assez voisin de 1, et si T est assez grand, cette 
probabilite est 0{N~^') quand A^ — t- 00 pour tout 6' < 5. 
En conclusion, si Ton choisit 

(7') |0| ^4(l + 5)(r2 + l)logAr, 6>0, 

si T est assez grand, la probabilite de (jrvl) 1 — 0{N~^ ) (A^ — t- 00) pour 
tout 5' < 5. II en est de meme pour la probabilite d 'avoir (jl]) pour toutes les 
fonctions x{t) portees par T points de G. 

C'est la variante aleatoire que j' avals annoncee. 

Le theoreme de Candes, Romberg et Tao dit que, etant donne une fonction 
x{t) portee par T points de G, si Ton choisit 

(8') |0| ^ 22(1 + (5)r log A^, 5>0, 

la probabilite d'avoir (jl]) pour cette fonction est 1 — 0{N~^) (voir [Can06| 
Th. 2.1] ; d'apres jCRTOGI (3.34)] on s'attendrait a 23 au lieu de 22 dans Q). 
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Dans leurs enonces, les auteurs n'insistent pas sur le fait que I'evenement dont 
ils estiment la probabilite depend du signal x{t), mais il en est bien ainsi dans 
la demonstration. 

L'interet de la version deterministe, outre la facilite de sa demonstration, est 
de montrer la portee possible en theorie du signal du probleme de la concen- 
tration des idempotents dans Z^r. 

La condition (17^ de la version aleatoire est plus exigeante que ([s]) des que 
T ^ 6. L'interet de la version aleatoire est qu'une fois il. choisi, (0) est valable 
pour tons les signaux portes par T points ; I'evenement dont on estime la 
probabilite ne depend que de T. Le prix a payer est de faire intervenir au 
lieu de T. 

Pourrait-on remplacer par une quantite o{T'^) (T — )■ oo) ? Si I'on s'at- 
tache a la probabilite de (jrvl) . la reponse est negative. En effet, ()iv|) entraine 

N-l\ ,^,iV + 4r2-l 



N\n\ = J2\Kit)\'^KiO)\l + ^) 

done 



42^2 J I I 42^2 

tGG 



9 N 



+ 4r2 - 1 ■ 

Peut-on retrouver et etendre le theoreme CRT en adaptant la demonstration 
de la version aleatoire ? La reponse est positive, mais il serait trop long de la 
donner ici. 

Exercices 

(a) Montrer que, pour la valeur choisie pour a, 

E -fc! E ( cos(27rnt/iV ^) < (t G \ {0}, ^ reel), 

done 

— tN——: ) (a = cosvr/i^). 

4T^ + 2 / 

(b) Quand T = 1, la donnee de x(l) determine x. Montrer que cependant 
on ne pent pas supprimer log ou le remplacer par une fonction de moins 
vite croissante dans la for mule (jl^), resp. ([l^) (prendre pair). 

(c) Expliciter la probabilite d'avoir (j!]) par le procede indique quand T = 2. 
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